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ABSTRACT 

Recent imaging campaigns indicate the likely existence of massive planets (~ 1 — 10 
Mj) on ^ 1000 year orbits about a few percent of stars. Such objects are not easily 
explained in most current planet formation models. In this Letter we use ensembles of 
100 N-body simulations to evaluate the potential for planet scattering during relaxation 
of dynamically active systems to explain the population of giant planets with projected 
separations up to a few 100 AU. We find that such a mechanism could indeed be at play, 
and that statistical samples of long period planets could place interesting constraints on 
early stage planet formation scenarios. Results from direct imaging and microlensing 
surveys are complementary probes of this dynamical relaxation process. 

Subject headings: planets and satellites: formation — methods: N-body simulations 



1. Introduction 



The ongoing discovery of exoplanets is providing a wealth of opportunities to test models of 
planetary syste m formation, and has already revealed that our own solar system is not necessar- 
ily typical (see lUdrv Santoa (|2007l ) for a review). Due to the natural biases of radial velocity 
and transit detection methods a major focus has been on short period planets, of all ages. Direct 
imaging efforts ar e biased oppositely , towards the detection of long period, young, and massive 
planets (see, e.g., iBeuzit et al.l 120071 . for a review). Microl ensing search es are also sensitive to 
longer period objects, independent of system age (see, e.g., iGaudil 120081 . for a review). Recent 
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direct imaging efforts have yielded constrai n ts on the likelihood of certain class e s of long pe 



riod gas-giant planets (jMasciadri et al.l l2005l : iBiller et al.l 120071 : iKasper et alj 120071 : iNielsen et al 



20081 ) . together with a number of increasingly robust claims of detection of l ong period planets 



2008 



or sub-stellar companions (e.g.. iNeuhauser et al. l2005l : iMohantv et al.l 120071: iQppenheimer et al 



Marois et al 



Lafreniere. Jayawardhana. &: van Kerkwiik 20081 ): and in particular 



Kalasetal 



(12003) and 



(120081') ■ These objects have projected separ ations of ~ 20 to a few 100 AU. The 



result of lLafreniere. Jayawardhana. &: van Kerkwijkl (120081 ) (LJvK08) is interesting, since it claims 
detection of an S'^'l Mj object at apparent separation of 330 AU from the parent star (IRXS 
J160929. 1-210524) - an approximately solar-mass star in the young (~ 5 Myr) Upper Scorpius 
association. Taken at face value this represents one positive detection out of ~ 80 potential planet 
hosting stars, or a detection rate of ~ 1% - although clearly this is a very crude evaluation since 
the detectability of companions is a complex function varying from star to star. 

Most protostellar/proto-planetary disks are approximately 100 AU in radius, and core ac- 
cretion models predict giant planet formation close to the water-ice snowline (< 10 AU), where 
coalescence timescales are close to, o r within, the disk lifetime against photoevap oration (e.g. 
Hollenbach. Yorke. Johnstond (j2000l )). Gravitational instability models (|Bossl 119971 ) might have 
a higher probability of forming massive planets at large radii - but this would require quite unusu- 
ally massive and extended disks. As pointed out by LJvK08, this suggests that migration and/or 
scattering mechanisms might be necessary to explain a planet at projected separation of 330 AU. 
Another possibility is that highly mass-imbalanced wide binary systems could form, with stellar 
and sub-stellar companions. 

In this Letter we present the results of numerical simulations of ensembles of 100 plane- 
tary systems that capture the effect of strong dynamical evolution through planet-planet scat- 
tering. We present a number of proof-of-principle predictions for objects with large semi-major 
orbital axes (> 100 AU), in particular those in young systems. Our methodology stems primar- 
ily from the work oflJuric fc Tremaind (1200811. and oth e r investigations of dynamical effects ^ (e.g. , 
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Juric &: Tremaind ([200a) suggest that, for the pur- 
poses of studying dynamical evolution, planet formation may be simplified into two stages. The 
initial episode (Stage 1) involves the assembly and evolution of a circumstellar and proto-planetary 
disk over ~ 10^ — lO'^ years. During Stage 1, planetary embryos are formed, destroyed, and inter- 
act strongly with their disk surroundings. Stage 2 takes place once the major part of the gaseous 
proto-planetary disk has dissipated. It involves the dynamical evolution of any planets that remain, 
including collisions, mergers, scattering, and significant planet ejection from the systems over the 
following tens to hundreds of million s of years. Using l a rge en sembles of statistically identical ini- 
tial conditions for planetary systems I Juric &: Tremaind (j2008l ) demonstrated that, in quite general 
circumstances, aspects of the observed distribution of giant exoplanet orbital eccentricities can be 
reproduced. This result is not strongly dependent on initial conditions as long as systems are 
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dynamically active, or partially active. It also demonstrates the utility of studying the statistical 
dynamical properties of a population of planetary systems, rather than trying to identify a physical 
process that prescribes orbital characteristics in any single system. 

Planets in dynamically active systems experience strong mutual interactions and frequent 
encounters and undergo relaxation. T his typically occurs when the spacing relative to the mutual 
Hill Radii (e.g., iHenon &: PetitI Il986l ) is small. Partially active systems contain only a subset of 
objects undergoing strong and frequent interactions, and while inactive systems still evolve, there 
are typically few or no major encounters in this case. Both active and pa rtially active systems ten d 
to lose planets as a function of time (instability timescales are given bv IChatteriee et al.l (|2008l )). 
This occurs through collisions (and mergers), capture by the parent star, or very often through 
complete ejection from the system. Thus, the Stage 2 dynamical paradigm predicts both a large 
population of ejected planets and a systematic excess of planets in very young systems compared 
to older systems. 



2. Simulations and methodology 



We ernploy the publicly available N-body integration routines in the MERCURY6 package 
(jChamberslll999l ) to follow the dynamical evolution of statistically identical ensembles of planetary 
systems. MERCURY6 also allows for close encounters and fully inelastic (no fragmentation) star- 
planet or planet-planet collisional mergers. We have found that a high-accuracy Burlisch-Stoer 
integrator typically provides relative energy errors of less than 10~^ up to integration times of 10'' 
years for the systems we are modeling. For simplicity we therefore use the Burlisch-Stoer BS2 
integrator in MERCURY6 and restrict our integrations to 10^ years - commensurate with the ages 
of the young plan etary systems we hope to place constraints upon. Our methodology is otherwise 
similar to that of lJuric &: Tremaind (120081 ). We note that for dynamically active systems all but 
the final ~ 10% of planet ejection/loss events have occurred by 10^ years. Our simulations do not 
include the effects of tidal dissipation between planetary objects and the parent star. Dissipation 
due to star-planet tides can act to dampen orbital eccentricities, and is an important physical 
process for planets with short periods. We choose to circumvent these issues by ignoring planets 
with semi-major axes of less than 0.1 AU in our final results. This is well motivated in that tidal 
damping is expected to operate on timescales ranging from ~ 10^ to ~ 10^2 

wars for objects 

with ~ 20 day orbital periods around solar mass stars (e.g., iMatsumura et al.ll2008l . and references 
therein). 

We adopt a (somewhat arbitrary) physical ejection radius of 2000 AU, but note that for semi- 
major axes beyond^ 1000 AU, secular effects from Galactic tides may become important (see, e.g., 
Higuchi et al.l 120071 ). It is possible that these tides would cause some of the planets we label as 
ejected to remain on bound orbits at very large separations. We neglect this possibility altogether 
in the present study. 
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2.1. Initial populations of objects 



We perform two sets of runs, each with statistical ensembles of 100 systems. We populate 
the systems with 10 planets and adopt initial properties such that the systems are dynamically 
active. Masses are drawn randomly from a uniform distribution in log(M) and are limited to a 
range of 0.1-10 Mj. Initial orbital inclinations (i) for all runs are drawn randomly from a Rayleigh 
(Schwarschild) distribution - where the peak, or mode, is equal to the distribution width parameter 
a, and we set cTj = 0.3 degrees. Orbital eccentricities are also drawn from a Rayleigh distribution 
with (Te = 0.3. In one run, that we denote as i?3o, the initial semi-major axes are drawn randomly 
from a uniform distribution in log{distance) between 0.1 AU and 30 AU. In the second run, denoted 
as i?ioO) the outer truncation is increased to 100 AU. All systems contain a central stellar mass of 
1M0. 

The i?3o and iiioo runs illustrate as simply as possible the dependence on initial conditions of 
the dynamical relaxation process. The two leading scenarios for planet formation, core accretion 
and direct gravitati onal collapse, should lead to significantly different initi al orbital outcomes for 



young planets (e.g., Lissauei 1993 : Boss 1997 JPaDalo izou &: Tera 



uem 



_ , 20061') . Each of these scenar- 
ios faces challenges (e.g., iMatzner Sz Levin! 12005. : .Rafikovll2005l : Dong et al. 2008 ) and they could 
in pr inciple both operate in individ ual systems or in the proto-planetary population as a whole 
(e.g.. lRibas &: Miralda-Escudd 120071 ). Given these complications, what is important in the present 
context is that the late dynamical architecture of exoplanetary systems could potentially allow us 
to discriminate between competing formation scenarios. Specifically, formation via core accretion 
is favored in the vicinity of the snow-line in the proto-planetary disk, typically well within a few 
tens of AU, while gravitational collapse is favored at larger orbital distances. This is the primary 
motivation for the two types of runs R^q and .Rioo performed here. 



2.2. Monte Carlo "snapshots" 

In order to model apparent separations of planet-star systems from our simulation outputs 
we must take into account the orbital characteristics and potential orientation of systems relative 
to a hypothetical observer. We adopt a simple algorithmic approach. First, the semi-major axis, 
eccentricity, and orbital period are used to construct a probability distribution for the instantaneous 
radial star-planet distance, assuming an observation made at a random time during the orbital 
period. We then assume that the orbital plane is randomly orientated with respect to the observer 
and project the star-planet radial separation onto the sky to obtain an apparent separation. 

We can draw an arbitrary number of "snapshot" projections for all remaining bound planets in 
our simulations. Since this step is not computationally intensive, the actual number of snapshots is 
simply determined by the need to adequately populate the mass - projected separation plane. Here 
we generate 1000 such projections per individual planet, irrespective of any other bound planet in 
the system. For simplicity, we then generate a cumulative distribution of projected separations for 
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all individual planets, ignoring entirely conditional probabilities associated with multiple planetary 
systems. This is sufficient for our present purpose, focused on the apparently most distant planet, 
but a more careful analysis would clearly require conditional probabilities for projected separations 
to be accounted for. 

Uncertainties in the separation distributions are estimated from the dominant Poisson statistics 
of the original simulations with a total of only 10 x 100 statistically-identical planets, as opposed 
to the more numerous projection realizations. For instance, if uniquely tagged planets, out of 
the total 1000 original ones, contribute to the statistics beyond some fixed projected separation, an 
uncertainty ^/N /N is applied to the apparent separation distribution. Plotted error bars/envelopes 
are la. 



3. Results 

In Figure 1 we plot the mean number of planets per system as a function of time for the two 
runs iiso and i2ioo- The more closely packed planets in i?3o result in a faster decay, and by 10^ 
years, when both runs are much more dynamically relaxed, the mean number of remaining planets 
in i?3o is slightly lower than that in R^nn, with ^ 2 planets pe r system in both cases. These results 



are in good agreement with those of lJuric &: Tremaind (|2008l ). who demonstrated that integrating 



to 10* years (by employing a symplectic integrator) results in only moderate further evolution, 
owing to the loss already of an average 80% of the initial population by 10^ years. 

By 10^ years in run i?3o, 25% of the initial population have collided with the central star, 12% 
undergo collision and merger with other planets and 44% are ejected (the majority by scattering 
into an unbound e > 1 orbit, some by passing beyond 2000 AU, see below). For i?ioo 21% collide 
with the central star, 11% undergo collision and merger, and 44% are ejected. This indicates that 
the primary difference in final planet numbers at 10'' years between R^q and -Rioo is almost entirely 
due to the higher rate of collision with the central star in R^q . 

In Figure 2 the eccentricity versus semi-major axis is plotted for all objects at several timeslices 
for the i?ioo run to illustrate the orbital evolution of the population. In the first ~ 10^ years a 
significant number of objects are scattered onto high eccentricity and large semi-major axis orbits. 
The great majority (> 99.2%) of objects that reach eccentricities e > 0.9 become unbound (e > 1) 
due to scattering close to periastron. 

In Figure 3 we plot the cumulative number of planets beyond a given apparent separation. 
In both iJso and .Rioo runs, some outward diffusion of the planet population is clearly seen by 
comparison to the initial distribution. Cumulative distributions after 10^ and 10^ years are shown, 
and only modest evolution is seen at these late times. In both runs, there is some reduction in the 
fraction of objects between a few AU and ~ 100 AU, from 10^ to 10'' years. 



For comparison with the putative planet detection of LJvKOS, we find that the population of 
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planets with 330 AU apparent separation or greater in the -Rioo run is 0.9(ib0.6)% with a variation of 
only ~ 0.2% between 10^ and 10^ years. In the R^q run, the corresponding numbers are 0.6(ib0.3)% 
at 10^ years and 0.2(±0.2)% at 10'' years. We note that these results are broadly consistent with 



Veras. Crepp. &: FordI (|2008l ) who find that ~ 0.5% of dynamically evolved (scattered) systems at 



late times (> 10^ years) could still harbor a long-term stable planet with e > 0.8 and a > 1000 AU. 

For a minimum apparent separation of 100 AU, at 10^ years, these fractions increase to 
7.1(±2.5)% in the i?ioo runs and 2.3(±1.0)% in the R30 runs. Later, at 10^ years, there is lit- 
tle evolution in the i?ioo results, but a drop to 1.4(ib0.6)% in the R30 results. 

Although these numbers are not suitable for a di rect comparison to c urrent imaging survey 
constraints, we nonetheless note that the constraint of Nielsen et al. (j2008l ) that less than 20% of 



systems should harbor > 4Mj planets between 20 and 100 AU is broadly consistent with the results 
of our simulations. The R^q runs indicate that less than ~ 16% of planets of all masses would be 
expected at larger than 20 AU apparent separations (at 10'' years). The corresponding i?ioo result 
is less than ~ 30%. 

The results shown in Figure 3 do not discriminate between planets of different masses as a 
function of projected separation. This is a particularly important issue for direct imaging surveys, 
with yields depending on a combination of ages, masses and projected separations. We believe 
that a larger number of statistical realizations than performed here is needed to address this point 
reliably and postpone a more detailed investigation to future work. 



4. Discussion 



Dynamical relaxation of young planetary systems can produce a population of long period 
planets within 10^ — 10'' years of planet formation. 

The results presented here are subject to a number of limitations. Although the outcome of 
dynamical relaxation of active systems in terms of ecce ntricities is not strongly dependent on the 
initial conditions set during Stage 1 planet formation (jjuric &: Tremaind l2008l ) . the precise rela- 
tionship between other population characteristics at ~ 10^ years in Stage 2 and initial parameters 
needs further exploration. For example, the differences between our R^q and i?ioo runs indicate 
that information on semi-major axis distributions is indeed retained in the population of objects 
at later times. The population of long period planets could therefore potentially provide important 
constraints on the outcome of Stage 1 planet formation. 

Furthermore, to make meaningful comparisons with observational surveys, one must allow for 
at least two biases: 1) Not all stars surveyed will necessarily form a population of giant planets 
during Stage 1. 2) N ot all systems with gian t planets will be dynamically active at the start of 
Stage 2. For instance. I Juric &: Tremaind (|2008l ) estimate on the basis of the observed distribution of 
exoplanet eccentricities that perhaps 75% of currently known planetary systems have been dynami- 
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cally active. Less dynamically active systems, which were not directly addressed here, may provide 
better records of the planet formation outcome but at the expense of a smaller range of projected 
separations (with respect to the dynamically relaxed apparent separations shown in Fig. 3). 

The detailed strategy of an imaging survey is also important for making meaningful compar- 
isons with theoretical expectations. If the planet formatio n phase (Stage 1) lasts between 10^ an d 



10^ years, the timescale for disk dispersal/evaporation (e.g.. lHollenbach. Yorke. &: Johnstonell200d ) . 
then a "blind" survey of young stars - where the target objects span a wide range of ages up to a 
few 10^ years - may not be well suited to probe the earliest Stage 2 phase. Indeed, the fast (Stage 
2) relaxation phase (Figure 1), relative to a longer duration Stage 1 gaseous phase, suggests that 
it may be statistically difficult to identify the systematically larger numb er of planets exp e cted in 



younger, dynamically-active proto-planetary systems. On the other hand, iThommes et al.l (|2008al ) 
present evolutionary scenarios based on mean-motion resonance locking and remnant planetisimal 
disk evolution in which the transition from Stage 1 to the Stage 2 relaxation occur much more grad- 
ually. Finding a larger number of planets in younger systems would provide strong observational 
evidence in support of the dynamical relaxation scenario. 

We conclude that the detection of massive planets at apparent separations of a few 100 AU need 
not be a challenge to the core accretion model of planet formation if dynamical relaxation occurs 
through extensive planet-planet scattering. Furthermore, the number of very long period planets 
resulting from scattering during the first few Myrs following disk dispersal may be consistent with 
current upper limits, and claimed detections, from all existing planet imaging campaigns. 

These are encouraging but very preliminary results. They suggest that a more systematic 
exploration of dynamical relaxation outcomes for long-period planets, as a function of initial con- 
ditions and system age, could provide a valuable additional means to interpret the existing and 
future yields of direct imaging surveys. Furthermore, since microlensing surveys probe the same 
projected planetary separations, but typically at much later system ages, they provide very useful 
complimentary information on the late stages of the dynamical relaxation process. Together such 
data would allow consistency checks to be performed on the systematic trends expected between 
early and late stages of orbital evolution in this promising model. 



The referee is thanked for comments that helped improve the paper. CS acknowledges the 
support of Columbia University's Initiatives in Science and Engineering, and NASA Astrobiology: 
Exobiology grant # NNG05GO79G. 
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Fig. 1. — The mean number of remaining planets per system is plotted as a function of time for 
runs i?3o and i?ioo- Most collisions (planet-planet and planet-star) and ejections occur within the 
first 10^ — 10^ simulated years. 



- 11 - 



1 

0.8 
0.6 



£4 O&o°„ VooSloft" 



I I iidiiii — m 




yr 

'rii<fiiiil I 



1 

0.8 
0.6 



10 100 1000 

a(AU) 



I MllillU d„l 



o oo <(J0»„ o o 
,0 o O Oo^nS* 
O » SO go" 



0.4 - ^o: 

_oV ----- 

0.2 



oo<b 



<D O" OOS 

o« ooOg 




0.1 



'«!)ogi« 10 yr - 

lit Ilftl'll °loiqi°Hll UJ 



1 10 100 1000 
a(AU) 




° I I P I I r I Piiffil UJ 




0.1 



1 



10 100 1000 



1 

0.8 
0.6 
0.4 
0.2 



I I iiiiiii — I I iiiiiii — loi I mill 



oo*8 o„ o ^o#„ 




0.1 



0„ ° °0 tfP O'^Oo „ 

fT O o°0 O O 
030 ° 0° o oo-o 

tt loll I Mill I U. 



1 



10 100 1000 



1 

0.8 
0.6 
0.4 
0.2 




I I iiMiii — rr 



TTIl — I I \*\\L,„ 

o o 



ntmi 




>^o= ,| %S o 4 



°l I I l^lTll L_L 



5«L 



P fli UJ 



0.1 



1 



10 



100 1000 



1 

0.8 
0.6 
0.4 

—I 0.2 I— oi? .o= 







' o< 

O S» J o" O O 

o ■ o 
ooo o « * 

°o -^oV „: 

'000*8 



8. n» 



,0 o Sof„°o^.s;t> 10 



I 1 1 



All f 1 1 



_UJ 



0.1 



10 100 1000 



Fig. 2. — The orbital eccentricities and semi-major axes of objects plotted at a number of times 
for the -Rioo run. Starting at the upper left and moving right, plots show the distribution at 0, 10^, 
and 10 years. Lower panels, left to right, arc for 10^ 10^ and lO'^ years. Of the highest eccentricity 
objects (e > 0.9) at any given time, less than ~ 0.8% remain on bound orbits by 10^ years. 
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Projected separation (AU) Projected separation (AU) 

Fig. 3. — The percentage of the remaining planet population with apparent separation greater than 
a given value is shown for the two models explored. Left panel: i?ioo ensemble results. Black curve 
and green error bars correspond to the population after 10^ years, blue curve and grey error bars 
correspond to the population after lO"^ years. Errors are computed from the Poisson uncertainty in 
the underlying number of planets used to seed the Monte Carlo projected separation distributions. 
The red dashed curve shows the initial conditions. Right panel: R^q ensemble results. In both 
panels, the consequences of an outward diffusion of planets from their initial distribution are clearly 
seen. Planets in the tail of the projected separation distribution typically have eccentricities e ^ 0.4 
(see Figure 2). 



